

// Table A3 
use cfps_pref, clear
foreach var of varlist qm1004 qm110*  {
	reghdfe `var' c.Q2##c.female c.Q3##c.female c.Q4##c.female, absorb(ybirth provcd) cluster(prov_ybirth)
	test 	Q2 Q3 Q4
	test 	c.Q2#c.female c.Q3#c.female c.Q4#c.female
} 

// Table A8
use cfps_food, clear
merge m:1 month using "D:\zhc\research\Labor and Health\birthmonth\workdata\month_var_pca"
drop _merge
egen std_food = std(food)
reghdfe log_food_expd2 sq2 sq3  sq4  if urban==0   , absorb(provcd urban) cluster(countyid)
test sq2 sq3  sq4
reghdfe log_food_expd2 std_food   if urban==0   , absorb(provcd urban) cluster(countyid)
sum log_food_expd2 if urban==0

// Table A12
use cfps_birthweight, clear
sum birthweight if male==0
sum birthweight if male==1
reghdfe birthweight Q2 Q3 Q4   if male==0, absorb(provbirth ybirth) cluster(prov_year_birth)
reghdfe birthweight Q2 Q3 Q4   if male==1, absorb(provbirth ybirth) cluster(prov_year_birth)
foreach y of varlist birthweight {			
	qui reg `y'  Q2 Q3 Q4  i.provbirth i.ybirth  if male==1
					est store `y'_m
	qui reg `y'  Q2 Q3 Q4  i.provbirth i.ybirth   if male==0
					est store `y'_f
	quietly suest `y'_m `y'_f,  vce(cluster prov_year_birth) 
	test ([`y'_m_mean]Q2 = [`y'_f_mean]Q2) ([`y'_m_mean]Q3 = [`y'_f_mean]Q3) ([`y'_m_mean]Q4 = [`y'_f_mean]Q4)  				
}

